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ABSTRACT 


This  report  describes  an  in-situ  technique  to  estimate  the 
following  parameters  of  a  phased-array  antenna: 

1.  the  relative  array-element  excitation  voltages, 

2.  the  array-element  tuning  phases,  and 

3.  the  RF  phase  shifts  at  the  array  elements. 


This  technique  has  several  significant  features.  First,  it  involves  the 
use  of  two  auxiliary  antennas.  One  is  a  remote  CW  source  directed  at  the 
phased-array  antenna.  The  other  is  a  passive  antenna  mounted  close  to 
the  phased-array  antenna.  Its  output  is  used  to  produce  a  reference 
phase  for  phase  measurer, lent  s .  Second,  it  contains  a  technique  to  reduce 
the  errors  in  phase  estimates.  Third,  it  takes  note  that  beam  steering 
uses  phase  sums  of  the  form  (<t>k+  where  ^  is  the  tuning  phase  for 

the  k-th  array  element  and  is  an  RF  phase  shift  of  the  array  element, 

and  pays  special  attention  to  reduce  the  errors  associated  with  the  esti¬ 
mates  of  these  sums.  Fourth,  it  assumes  the  use  of  a  reasonably  stable 
and  strong  CW  source  of  commercially  available  quality.  No  other  assump¬ 
tions  are  made. 


Experimental 
phased-array  antenna 


results  obtained  with  a  295-element  S-band  space-fed 
are  given. 
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I.  INTRODUCTION 

A  phased-array  antenna  requires  periodic  measurements  of  the 
following  parameters: 

1.  the  relative  array-element  excitation  voltages, 

2.  the  array-element  tuning  phases,  and 

3.  the  RF  phase  shifts  at  the  array  elements. 

These  parameters  are  needed  to  compensate  for  known  and  unknown  varia¬ 
tions  in  the  array-element  characteristics,  the  antenna  geometry,  and  the 
antenna  orientation.  This  task  is  preferably  carried  out  in  the  back¬ 
ground,  so  that  it  does  not  interfere  with  the  normal  operations  of  the 
antenna.  The  measurement  technique  is  necessarily  in-situ,  and  should  be 
fast,  interruptable,  and  not  involve  the  use  of  an  aperture  probe  if  the 
antenna  is  to  be  returned  to  operation  in  small  fractions  of  a  second. 

Several  techniques  which  can  be  used  to  estimate  the  above  array- 
element  parameters  under  the  conditions  specified  have  been  published  in 
the  open  literature.  Notable  examples  are  those  described  in  the  papers 
by  Lowenschuss  [1],  Huschelrath  and  Sander  [2],  and  Alexander  and  Gray 
[3].  Lowenschuss  uses  an  auxiliary  antenna  in  the  far-field  zone  of  the 
array  antenna.  A  cable  links  this  source  with  the  array  receiver  to  pro¬ 
vide  a  reference  phase  for  phase  measurement.  Three  assumptions  are 
made:  (i)  the  nominal  180*  RF  phase  shift  of  each  array  element  is  close 

to  the  actual  value;  (ii)  the  change  in  excitation  voltage  at  an  array 
element  can  be  neglected  when  a  phase  bit  is  changed;  and  (iii)  the  tuned 
-tate  of  an  array  element  can  be  identified  as  one  of  the  phase  states  of 
the  array  element.  The  RF  phase  shift  estimation  technique  in  [1]  has 
been  modified  for  estimation  with  antennas  which  contain  thousands  of 
elements  by  Kahrilas  and  Jahn  [4]  and  by  Blake,  Schwartzman,  and  Esposito 
[5].  The  descriptions  of  these  modified  techniques  are  very  sketchy. 
Huschelrath  and  Sander  [2]  use  essentially  the  same  reference  phase 
generation  technique  as  in  [1].  They  assume  very  accurate  knowledge  of 
the  distances  between  the  auxiliary  antenna  and  the  individual  array  ele¬ 
ments.  The  technique  is  not  described  in  detail  in  the  reference. 
Alexander  and  Gray  [3]  use  an  auxiliary  antenna  in  the  far-field  zone  cf 
the  array  antenna.  They  do  not  use  a  cable  linkage  between  the  auxiliary 
antenna  and  the  array  radar.  Instead,  when  an  array  element  is  tested, 
the  reference  phase  is  identified  as  the  phase  of  the  sum  contribution  of 
all  the  other  array  elements.  The  three  assumptions  in  [1]  are  also 
made. 


There  are  both  favourable  and  unfavourable  features  in  each  of 
the  above  reference  phase  generation  techniques.  The  techniques  in  fl] 
and  [2]  produce  reference  phases  which  are  almost  perfectly  synchronized 
with  that  of  the  signal  at  the  auxiliary  antenna.  This  is  favourable  for 
RF  phase  shift  measurements.  However,  the  auxiliary  antenna  in  these 
techniques  cannot  be  located  far  away  from  the  array  antenna.  This  is 
unfavourable  for  tuning  phase  measurements,  because  accurate  distance 
measurements  are  needed  to  correct  for  wavefront  sphericity.  The  tech¬ 
nique  in  [3]  uses  an  auxiliary  antenna  which  can  be  located  at  a  large 


distance  away  from  the  array  antenna.  This  is  favourable  for  tuning 
phase  measurements.  However,  the  reference  phase  in  this  technique  is 
derived  from  the  antenna  output  signal.  This  is  unfavourable  for  RF 
phase  shift  measurement,  because  any  error  in  reference  phase  is  added  to 
the  RF  phase  shift  estimate. 

The  three  assumptions  made  in  F 1 ]  and  f 3 ]  are  usually  not  satis¬ 
fied.  It  is  not  uncommon  to  have  a  difference  of  more  than  a  few  degrees 
between  the  nominal  180°  RF  phase  shift  and  the  actual  value.  It  is  also 
not  uncommon  to  have  a  change  of  more  than  10%  in  the  excitation  voltage 
at  an  array  element  when  its  phase  state  is  changed.  A  10%  change  in  ex¬ 
citation  voltage  alone  can  lead  to  msny  degrees  of  error  in  RF  phase 
shift  estimates.  The  third  assumption  on  tuning  phase  is  valid  only  if 
the  number  of  phase  bits  in  the  phase  shifters  is  large.  In  the  special 
case  of  3-bit  phase  shifters,  for  example,  the  tuned  state  of  an  array 
element  can  deviate  from  the  phase  state  of  the  element  by  up  to  22.5* 
(=180°/2  The  technique  in  [2]  is  difficult  to  evaluate  and  implement, 
because  not  enough  details  are  published  in  the  paper. 

This  report  describes  a  new  technique  to  measure  the  relative  ex¬ 
citation  voltages,  tuning  phases,  and  RF  phase  shifts  of  the  elements  of 
a  phased-array  antenna.  It  has  several  significant  features.  First,  it 
incorporates  a  new  reference  phase  generation  technique  which  has  all  the 
favourable  features  in  [1],  [2],  and  [3]  and  none  of  their  unfavourable 
features.  Second,  it  takes  note  that  the  phase  states  of  an  array  ele¬ 
ment  are  related  to  each  other  and  uses  this  property  to  reduce  the  er¬ 
rors  in  the  phase  estimates.  Third,  it  takes  note  that  beamforming  uses 
phase  sums  of  the  form  (^k^kA)*  where  ^  is  the  tuning  phase  for 
the  k-th  array  element  and  is  an  R*'  phase  shift  of  the  array  ele¬ 
ment,  and  pays  special  attention  to  reduce  the  errors  associated  with  the 
estimates  of  these  sums.  Fourth,  it  assumes  the  use  of  a  reasonably 
stable  and  strong  CW  source  of  commercially  availble  quality.  No  other 
assumptions  are  made. 

There  are  seven  sections  in  this  paper.  Section  II  describes  the 
array  antenna  used  to  develop  the  new  technique  and  the  relationship 
among  the  RF  phase  shifts  of  its  phase  control  modules.  Section  III  con¬ 
tains  the  derivation  of  the  equations  used  to  estimate  the  excitation 
voltages,  tuning  phases,  and  RF  phase  shifts.  Section  IV  discusses  the 
errors  associated  with  the  phase  estimates.  Section  V  describes  the  ex¬ 
perimental  setup  to  generate  a  reference  signal  from  a  remote  CW  source 
and  the  measurement  of  phased-array  antenna  output  signals.  Section  VI 
contains  two  examples  on  the  estimation  of  excitation  voltages,  tuning 
phases,  and  RF  phase  shifts.  The  first  uses  a  distortion-free  space-fed 
phased-array  antenna.  The  second  uses  a  distorted  antenna.  Section  VII 
contains  concluding  remarks. 

II  THE  PHASED  ARRAY-ANTENNA 

The  phased-array  antenna  used  in  the  experimental  verification  of 
the  estimation  technique  is  on  loan  to  the  Canadian  Department  of 
National  Defence  from  the  Royal  Signal  and  Radar  Establishment,  Malvern, 
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United  Kingdom,  under  the  auspices  of  the  The  Technical  Cooperation  Pro¬ 
gram.  It  is  a  reciprocal  S-band  antenna,  consisting  of  a  rectangular 
feed  horn  mounted  1.11  m  behind  a  vertical  planar  array  as  shown  in 
Figs.  1  and  2.  The  planar  array  contains  295  elements  configured  into  an 
eliptical  aperture  with  a  1.67  m  horizontal  axis  and  a  1.05  m  vertical 
axis.  In  each  array  element  is  a  3-bit  PIN-diode  phase-control  module 
attached  to  flared  waveguide  horns  at  the  ends.  The  PIN  diodes  are  used 
to  switch  three  transmission  line  segments  with  nominal  lengths  0.125, 


and 


hs  in 

and 

out  of  the 

signal 

path 

,  thereby  produc- 

ft  s . 

The 

RF  phase  shifts  at 

the 

k-th  array  ele- 

{ek£i 

*=0 , 

,1,2,.  ..,7}, 

.  They 

sat  i: 

sfy  the  relations 

6k*  = 

0 

» 

*=0 

(1) 

0k  *  ~ 

*  X 

45°  , 

*=1,2,4 

9 

(2) 

ek*  = 

8kl 

+  ek2  » 

*=3  , 

(3) 

8k*  = 

ek4 

+  8k(*-4)  » 

*=5,6,7 

. 

(4) 

The  independent  RF  phase  shifts  are  3ki,  3k2>  an^  &k4-  Their  actual 
values  depend  on  the  actual  lengths  of  the  nominal  0.125,  0.25,  and  0.5 
wavelength  transmission  lines,  respectively,  in  the  phase  control  module 
of  the  k-th  array  element. 

Ill  EQUATIONS 

In  this  section,  an  expression  for  the  antenna  output  due  to  a 
remote  CW  test  source  is  given.  The  conditions  for  the  measurement  of 
antenna  outputs  are  stated.  Equations  for  the  estimation  of  array- 
element  excitation  voltages,  tuning  phases,  and  RF  phase  shifts  are  then 
derived.  A  summary  is  given  at  the  end. 

A.  Expression  for  Antenna  Output  Amplitude 

The  complex  amplitude  at  the  phased-array  antenna  output  due  to  a 
remote  CW  source  can  be  written  as  a  sum  of  the  contribution  from  a 
single  specified  array  element  plus  a  second  term  comprising  all  the 
other  contributions, 


Ak*  *  ak*  ej(°k+Pk£)  +  Bk  .  (5) 

Here,  k  signifies  the  quantities  associated  with  the  specified  element 
and  *  signifies  the  phase  state  of  this  element.  The  elemental  amplitude 
from  this  element  has  an  effective  magnitude  ak£.  Its  phase  is  ak 
before  phase  shifting.  Amplitude  Bk  represents  the  sum  of  all  other 
contributions  to  Ak£.  It  includes  the  sum  of  elemental  amplitudes  from 
all  other  array  elements,  mutual  coupling,  and  background  noises. 


a30WnN  A\Od 


Fig.  1  Layout  of  array  elements  as  seen  from  the  back.  The  array  has 

295  elements  configured  into  an  elliptical  aperture  with  a  1.67  m 
horizontal  axis  and  a  1.05  m  vertical  axis. 


JMN  NUMBER 
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Fig.  2  Rear  view  of  the  space-fed  phased-array  antenna  with  the  feed 
horn  in  position. 
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shown  in  Figure  3.  Also  shown  in  this  figure  are  the  tips  of  other 
amplitudes  in  the  set  {a^;  4=1, 2,..., 7}  in  the  special  case  of  ak£  = 
ak  and  Bk£  =  4x45°.  These  tips  lie  on  a  circle  of  radius  ak  .  The 
complex  amplitudes  representing  the  contributions  from  the  k-th  array 
element  are  denoted  by  vectors  radiating  from  the  tip  of  Bk  and 
terminating  at  the  circumference  of  the  circle. 

Because  the  sum  of  ctk  and  Bkjj,  appears  in  the  expression  for 


Ak£  it  is  convenient  to  denote  this  sum  as 

Tk4  *  ak  +  6k4  • 

From  the  properties  of  the  in  (1)  to  (4),  one  has 

YkO  "  “k  >  (7) 

^k4  ~  Yko  *  &k4«  4=0, 1,2 . 7  ,  (8) 

Yk4  “  Tk(i-4)  *  ^k4 »  A-A.5,6,7  ,  (9) 


and 


Ykl  +  Yk2  *  Yk0  +  Yk3 


(10) 


B.  Assumptions 

In  the  rest  of  this  paper,  it  shall  be  assumed  that  ak£,  ak 
and  Bk  are  constants.  This  means  that  the  experimental  setup  must  be 
such  that  the  following  conditions  are  satisfied  during  amplitude 
measurement : 


1.  The  output  power  of  the  CW  source  is  sufficiently  stable  over 
the  duration  of  an  antenna  calibration  run. 

2.  The  phase  coherency  between  the  CW  source  and  the  reference 
signal  at  the  radar  receiver  is  very  high. 

3.  The  magnitude  of  fluctuations  in  the  noise  background  is 
small  compared  with  akjj,. 

The  first  condition  is  not  difficult  to  satisfy.  In  practice  a 
CW  source  with  an  output  power  which  varies  by  less  than  2X  over  the 
duration  of  a  calibration  run  is  good  enough  for  amplitude  measurements. 
This  type  of  CW  source  is  readily  available  commercially.  The  detailed 
experimental  setup  and  the  suppression  of  noises  to  satisfy  the  remaining 
two  conditions  are  described  in  Section  V. 

Amplitude  B^  includes  the  contribution  due  to  mutual  coupling 
among  the  array  elements.  Strictly  speaking,  it  is  also  dependent  on  the 
phase-shift  state  4  of  the  k-th  array  element.  However,  this  dependence 
on  4  is  very  weak  and  can  be  ignored  completely. 


REAL  AXIS 


Geometric  representation  of  in  the  Argand  diagram  and  the 
tips  of  other  complex  amplitudes  in  the  set  {a^j;  £«0, 1 ,2 } 
in  the  special  case  of  a^£  ■  constant  and  $kl  ■  i.x45*.  The 
tips  of  the  complex  amplitudes  in  the  set  lie  on  a  circle  centred 
at  the  tip  of  The  vector  drawn  from  the  tip  of  to  the 

tip  of  Afc$  is  the  elemental  amplitude  from  the  k-th  array 
element.  The  angle  between  this  vector  and  the  real  axis  is 
Tk*.  Yk*mak+ekt* 
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C.  Estimation  of  Relative  Excitation  Voltages 

The  value  of  ak£  in  (5)  is  dependent  on  both  k  and  £.  It  is 
dependent  on  k  because  the  coupling  between  the  feed  horn  and  an  array 
element  is  dependent  on  the  position  of  the  array  element.  It  is 
dependent  on  £,  because  the  insertion  loss  at  the  k-th  array  element  is 

dependent  on  the  signal  path  specified  by  £.  Only  the  average  value  of 

the  ajc£S  over  all  the  phase  states  of  the  k-th  array  element  is  of 
interest  here.  This  average  is  dependent  on  how  it  is  calculated.  One 

method,  suggested  by  Fig.  3,  is  to  calculate  it  as 

a 

ak  =  radius  of  the  best-fit  circle  through  the  set  of 

amplitudes  {a^;  £=0, 1 ,2 , . . . , 7  }  in  the  Argand  diagram,  (11) 

This  radius  is  also  the  radius  of  the  best-fit  circle  through  the  set  of 
elemental  amplitudes  |ajc£e^*c^;  £=0,1,2 . 7}. 

The  estimate  of  relative  excitation  voltage  at  the  k-th  array 
element  is  defined  here  as 


a 


k=l  ,2,3 . 295 


(12) 


where 


(13) 


is  a  normalization  constant,  K'  is  the  number  of  functioning  elements, 

and  l*  is  a  summation  over  the  functioning  elements.  Thus  the  bks 
defined  here  satisfy  the  condition 


!’bk  -  K'  , 

so  that  the  average  of  the  relative  excitation  voltages  over  all  the 
functioning  elements  is  unity.  An  array  element  is  arbitrarily 
identified  as  a  functioning  element  if  ak  is  larger  than  30Z  of  its 
theoretical  value.  This  choice  of  threshold  is  based  on  the  Section  VI 
results  to  be  presented  later. 

The  position  and  the  radius  of  the  best-fit  circle  through 

lAk£5  £■0,1,2 . 7}  are  dependent  on  the  criterion  foi  the  construction 

of  this  circle.  It  is  convenient  to  define  this  circle  as  the  circle 
which  minimizes  the  cost  function 

s  =  l  (dk£  *  ak)2 
£-0 


(15) 
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where  dk£  is  the  distance  of  Ak£  from  the  circle  centre  Bk,  an 
estimate  of  Bk.  This  cost  function  is  chosen  for  ease  of  analysis  and 
the  fact  that  it  leads  to  tractable  closed  form  solutions.  In  Appendix 
A,  it  is  shown  that  the  coordinates  of  Bk  can  be  calculated  in  closed 
form  as 

Re{Bk}  =  [f(x ,y)-g(x,y) )/h(x,y)  ,  (16) 

A 

Im{Bk}  =  [f(y ,x)-g(y,x)]/h(x,y)  ,  (17) 

where 


f(x,y)  =  ?.(y2-yy)(x3-xx2+xy2-xy2)  ,  (18) 


g(x,y)  =  2(xy-xy)(x2y-x2y+y3-yy2)  ,  (19) 

h(x,y)  =  4[x2-xx)(y2-yy)-(xy-xy)2]  ,  (20) 

and 

-  !  7 

xUyV  »  —  l  [Re(Ak£)]U[lm(Ak£)]V  .  (21) 

8  Jt=0 

The  radius  of  the  circle  is  given  by 

*  —j  —j  -  *  —  *  _  ^ 

ak  -  (x^+y  -2xRe(Bk)-2yIm(Bk)  + | Bk | ^ )  .  (22) 

D.  Estimation  of  RF  Phase  Shifts 


The  RF  phase  shifts  at  the  k-th  array  element  are  8ko»  8kl * 
8k2>***»8k7*  *n  t*ie  following,  these  phase  shifts  are  estimated 
through  the  estimation  of  the  Yk£S  given  by  (6).  There  are  five  steps 
in  the  procedure.  First  a  set  of  initial  estimates  denoted  by 


bio  ;  ^*o ,1,2 . 7}  is  constructed . 


Second,  the  estimate  of  Bk^  is 


calculated 
est imates 


.a  Third,  condition  (9)  is  imposed  to  produce  a  set  of  improved 
{Yki  »  4=0 ,1,2,. ..,7].  Fourth,  condition  (10)  is  imposed  to 


produce  the  final  estimates  bk£  5  4=0,1 ,2,. .. ,7}.  Finally,  estimates  of 
8ki  and  are  calculated.  Imposition  of  conditions  (9)  and  (10)  on 


the  estimates  of  the  Yk£S  leads  to  reductions  in  the  errors  in  phase 
estimates.  A  discussion  on  this  error  reduction  is  given  in  Section  IV. 


To  derive  an  initial  estimate  of  the  Yk£S,  it  is  noted  that  the 
phase  of  (Ak£  -  Bk)  can  be  identified  as  the  initial  estimate  of 
Yk4»  i *e • » 
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Y^A}=  phase  (Ak£-  ^1,  £=0,1 ,2, . . .  ,7  .  (23) 

There  are  four  values  of  phase  difference  given  by 

[ Yk ( £+4 )~  Yk£  I»  £=0,1,2, 3.  According  to  (9),  each  can  be  used  to  obtain 

an  estimate  of  8k4«  Because  the  nominal  value  of  8k4  is  180®,  these  four 
initial  estimates  of  8^4,  identified  with  £=0  to  3,  can  be  calculated  as 


s(1)m  -  C(1) 

Sk4  (A)  \(  £+4) 


-  Y(1)+ 
Tk£ 


nk£  x 


360® 


(24) 


where  njt£  is  an  integer  which  restricts  ^k^(£)  to  the  range  [0®,360") 

centred  at  the  nominal  value  of  8k4*  The  final  estimate  of  8^4  is 
defined  as  the  average  of  the  initial  estimates, 


Sk4 


1 

7T 


I  &><«> 

£=0 


(25) 


Condition  (9)  is  then  imposed  through  the  construction  of  a  new 

/  A  \ 

set  of  estimates,  denoted  by  ;  £=0,1 ,2 , . . . ,7}  defined  as 

A  A  A 

Ykl)+  0.5[4i)(£)-8U)],  £=0,1,2, 3  , 


v(2) 

Tk£ 


(26) 


Y^-  0.5[8^)(£-4)-8^)],  £=4,5,6, 7  . 

It  is  a  simple  matter  to  verify  that  these  improved  estimates  satisfy  the 
relation 


(2)  (2) 

Vi  -  k( £-4)*  ek4»  *“4’5’6’7 


(27) 


Relation  (10)  is  imposed  through  the  use  of  a  two-step 
calculation  procedure  derived  in  Appendix  B.  Initially,  a  dummy  phase 
angle  is  calculated  as 


6 


k 


(2) 

kl 


+  p^  x  360* 


(28) 


where  p^  is  an  integer  which  restricts  6^  to  the  range  [-180* ,180® ) . 

This  range  is  chosen  because  <5^  is  zero  if  (10)  is  satisfied.  The  second 

step  calculates  the  final  estimates  { Yk £»  £=0 , 1 ,2 , . . . ,7 }  as 
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A(2) 

Yl a  =  Yk£  + 


,1,2,..  . ,  7 


(29) 


The  e^s  are  derived  in  Appendix  B  and  are  given  by 


e4  = 


0.25, 

-0.25, 


4=0, 3, 4, 7, 
4= 1,2, 5, 6. 


(30) 


Estimates  of  RF  phase  shifts,  including  Bk0  anc*  &k4  »  are  calcu¬ 
lated  from  the  final  estimates  of  the  Yk£S  as 


8k£  =  Yk£  -  Yko  +qk£x360°,  4=0, 1 ,2 , . . . ,7  .  (31) 

A 

Here  qk£  is  an  integer  which  restricts  $k£  to  the  range  f (4x45°-180° ) , 
(4x45*+180°)) .  This  range  is  centred  at  the  nominal  value,  4x45°,  of 

0k4- 

A  A 

The  Bk4  calculated  with  (31)  is  the  same  as  the  Bk4  calculated 
with  (25).  The  0k£s  satisfy  relations  (1),  (3)  and  (4),  i.e., 


0k4  a  0  » 

4=0 

(32) 

AAA 

ek4  *  &kl  ♦  ^k2» 

4=3  , 

(33) 

0k4  *  Bk4  +  ^k(4-4), 

4=5, 6, 7 

(34) 

E.  Estimation  of  Tuning  Phase 

The  main  beam  of  the  array  antenna  is  pointed  at  the  remote  CW 
source  if  all  the  elemental  amplitudes  from  the  array  elements  arrive  at 
the  feed  horn  with  the  same  phase.  In  this  case,  the  radar  receiver 
output  signal  has  the  largest  possible  amplitude.  Angle  ak  in  (5)  is 
the  phase  of  the  k-th  elemental  signal  at  the  radar  receiver  output.  A 
phase  reduction  equal  to  a^,  equivalent  to  a  phase  advance  of 
(360*-^),  applied  to  the  k-th  elemental  signal  would  force  it  to  have 
a  zero  phase  at  the  radar  receiver  output.  If  all  the  elemental 
amplitudes  were  to  have  zero  phase  at  the  radar  receiver  output,  the 
radar  receiver  output  would  have  the  largest  possible  amplitude.  Hence, 
apart  from  a  constant  common  to  all  array  elements,  is  the  phase 
reduction  which  must  be  applied  to  the  k-th  elemental  amplitude  to  steer 
the  main  beam  of  the  antenna  at  the  remote  CW  source. 

T 

The  theoretical  value  of  ak ,  denoted  by  at^  ,  can  be  calculated 
from  the  antenna  geometry,  the  CW  source  direction  relative  to  the 
antenna  boresight  direction,  and  the  array-element  insertion  phases.  The 
difference 


is  the  tuning  phase  for  the  k-th  array  element.  From  (7),  an  estimate  of 
ak  is  given  by 

A  A 

ak  =  Yk0  .  (36) 

Therefore  an  estimate  of  4>k  is 

A  A 

♦k  =  Yk0  “  ak  •  <37) 

F.  Summary 

Let  {Ak£j  k=l ,2 ,3 , . . . ,295;  £=0,1 ,2 , . . . ,7 }  be  the  phased  array 
antenna  outputs  measured  with  a  remote  CW  source.  The  estimation  of 
relative  array-element  excitation  voltages,  tuning  phases,  and  RF  phase 
shifts  can  be  carried  out  in  eight  steps: 

Step  1  Set  k=0 

Step  2  Replace  k  by  k+1 

A 

Step  3  Calculate  circle  radius  ak  and  circle  centre  Bk  with  (16)  to  (22) 

1  \ 

Step  4  Calculate  {yk£  ;  £=0,1 ,2 ,. . . ,7}  with  (23) 

Calculate  {yk£  ;  £=0,1 ,2 , . . . ,7 }  with  (24)  to  (26) 

A 

Calculate  {yk£;  £=0,1 ,2 , . . . ,7  }  with  (28)  to  (30) 

Calculate  estimates  of  RF  phase  shifts  { 3k£ ;  £*0,1 ,2 , . . . ,7 } 
with  (31) 

Step  5  Calculate  the  estimate  of  tuning  phase  $k  with  (37) 

Step  6  If  k*295  proceed  directly  to  Step  8 
Step  7  Return  to  Step  2 

Step  8  Calculate  the  estimates  of  the  relative  array-element  excitation 
voltages  {bk;  k=l ,2 ,3, . . . ,295 }  with  (12)  and  (13). 

A  A  A 

The  estimates  ak,  ak,  and  the  0k£S  can  be  checked  in  two  places. 
First  the  8^^(£)s  calculated  in  (24)  must  be  close  to  180*.  Second,  the 


6^  calculated  with  (28)  must  be  close  to  zero,  because  its  expected 
value  is  zero. 


IV  ERRORS  IN  PHASE  ESTIMATES 

Let  the  initial  estimates  given  by  (23)  be  related  to  the  correct 
values  by 

Yk£^  =  \l  +  ek£  »  (38) 

k  =  1,2,3,... ,295  , 

£  *  0,1 ,2 . 7  , 

where  is  the  error  in  the  initial  estimate  of  Yk£.  The  statistical 
properties  of  the  e^s  have  been  studied  with  computer  simulations.  The 
results  show  that,  for  a  given  array  element,  the  errors  are  independent, 
have  zero  mean,  and  the  same  variance,  i.e. 


E{Gk£}  =  0  . 


0,  £*£  , 
EfcklEkt' 1  '  2 

k 


£,£'-0,1,2,. ..,7 


The  variance  of  ejt£  depends  on  k.  Usually,  it  is  smaller  if  the  element 
excitation  voltage  is  larger. 

From  (24)  to  (30),  one  can  show  that  the  estimates  of  the 
independent  RF  phase  shifts  are  related  to  the  correct  values  by 

A 

6kl  "  6kl  +  °-25^"ek0+ekl“ek2+ek3"ek4+ek5~ek6+ck7)  *  (41) 


8k2  *  &k2  +  °-25(-ek0~ekl+ek2+ek3"ek4"ck5+ek6+ek7>  >  <42> 


6k4  "  ®k4  +  °*25(”ek0"ekl~ek2"ek3+ek4+ek5+ek6+ek7)  •  <43> 

A  A 

These  expressions,  together  with  (39)  and  (40),  show  that  6^1,  6^2  and 
are  unbiased  estimates, 
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with  variances 


var^Jt}  =  0.5ok  £=1,2,4  .  (45) 

The  estimate  of  the  tuning  phase  is  related  to  the  correct  value 

<f>k  by 

A  A 

4>k  =  *k  +  ^k0_  ^k0> 

=  <j)k  +  0 . 25(  2eko+Gkl+ek2+ek4-ek7 )  • 

The  first  line  on  the  right  hand  side  of  this  expression  is  derived  with 
(7),  (35),  and  (37).  The  second  line  is  derived  with  (38)  and  (24)-(30). 
This  estimate  is  unbiased, 


a 

E{4>k}  -  4>k  »  (*7> 

and  has  variance 

varj^k}  =  0.5ak  .  (48) 

Phase  sums  {(4>k+8k£)»  k=l,2,3 . 295;  £=0,1  ,2, .  . .  ,7}  are  used  in 

beam  steering.  For  example,  to  steer  a  beam  at  the  remote  CW  source,  the 
phase  state,  identified  as  £(k),  at  the  k-th  array  element  must  be  such 
that  ( (J*k+&k£(k))  cl°9er  t0  ~aI  than  any  other  phase  sum  in  the  set 
{( V-&k£>;  *-0,1,2,. ..,7}.  The  estimates  of  the  phase  sums  can  be 
calculated  as 


A  A  A 

(41c  +  0k£>  -  41c  +  0k£  »  (49) 

k=l  ,2,3,... ,295 
£-0,1 ,2,. . . ,7  . 

Using  relations  (1)— (4),  ( 32 ) — ( 34 ) ,  and  (41 )— (46) ,  one  can  show  that  these 
estimates  are  unbiased. 


E{(4>k  +  &k£)  1  M  4>k  +  &k£ 


(50) 


and  have  variances 


var 


{(4k  +  &k£>} 


0.5oJ 


(51) 


For  example,  the  estimate  of  ($k+®k7)  can  be  calculated  with  (44), 


Experimental  setup  to  measure  the  output  of  the  phased-array 
antenna.  The  CW  source  was  slightly  over  one  kilometre  in  front 
of  the  phased-array  antenna.  The  dish  antenna  was  located  at  a 
distance  of  about  7  metres  from  the  phased-array  antenna.  Its 
output  was  used  as  a  reference  signal  for  the  quadrature  sampling 
of  the  phased-array  antenna  output. 


Step  1 

Set  up  the  experiment  as  shown  in  Fig.  5.  Steer  the  CW  source  at 
the  phased-array  antenna.  Steer  the  dish  antenna  at  the  source.  Steer 
the  phased-array  beam  away  from  the  source  to  make  the  initial  I  and  0 
channel  outputs  as  small  as  possible. 

St  ep  2 

Set  k=0. 


Step  3 

Replace  k  by  k+1. 


Step  4 

Cycle  the  k-th  array  element  through  all  its  eight  phase  states. 
Take  M  individual  antenna  output  amplitudes  at  each  phase  state.  Identify 
the  m-th  amplitude  at  the  *-th  phase  state  as  D^j^m). 

Step  5 

Return  the  k-th  array  element  to  its  original  phase  state,  i.e., 
the  phase  state  at  the  end  of  Step  1. 

Step  6 

If  k=295,  proceed  directly  to  Step  8. 

Step  7 

Return  to  Step  3. 

Step  8 

Calculate  Aj<£  as 

1  M 

Ak*  "  -  l  >  (56) 

M  m«l 


k*l ,2,3,... ,295  , 

*■0,1,2, . . . ,7 

A  choice  M*2000  was  made.  This  choice  produced  estimates  , 

A  A  . 

&k2»  8r4  ,  and  4>k  which  had  standard  deviations  less  than  0.25  degree  in 

repeated  measurements  with  the  central  element  in  the  array.  This  central 
element  was  identified  with  kM148  (see  also  Table  I). 

VI  RESULTS 

The  results  of  two  experiments  to  estimate  relative  excitation 
voltages,  tuning  phases,  and  RF  phase  shifts  are  presented  here.  The 
first  was  carried  out  with  a  distortion-free  phased-array  antenna.  The 
second  was  carried  out  with  the  feed  horn  rotated  to  provide  an  example  of 
calibration  with  a  distorted  antenna. 
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A.  Distortion-Free  Antenna 

Fig.  5  shows  the  best-fit  circle  through  the  set  of  complex 
amplitudes  jA^;  k=148;  £=0,1 ,2, . . . ,7 }  measured  with  the  central  element, 
identified  with  k  =  148,  cycled  through  all  its  RF  phase  states.  The  tips 
of  the  Aj^s  are  marked  with  crosses  labelled  with  £.  The  positions  of  the 
crosses  relative  to  the  circle  might  suggest  that  the  gain  for  the 
I-channel  output  (real  component  of  the  Aj^s)  should  be  increased.  This 
modification  of  the  I-channel  gain  was  not  carried  out,  because  the 
complex  amplitudes  measured  with  other  array  elements  did  not  support  this 
gain  modification. 

The  estimated  values  of  the  relative  excitation  voltages  are 
plotted  in  Fig.  6.  They  were  in  good  agreement  with  the  theoretical 
values.  The  voltages  in  the  centre  of  the  array  were  relatively  larger 
than  those  near  the  edges,  because  the  feed  horn  was  closest  to,  and 
pointed  at,  the  centre  element.  The  dots  in  the  figure  denote  the 
positions  of  the  defective  elements.  The  excitation  voltages  at  these 
elements  were  less  than  10%  of  the  theoretical  values.  The  excitation 
voltages  at  the  other  elements  were  more  than  50%  of  the  theoretical 
values . 


Listed  in  Table  I  are  some  RF  phase  shifts  obtained  both  by 
estimation  and  by  bench  measurements.  Array  elements  k=128,  148,  and  175 
were  in  the  centre  of  the  middle  column,  identified  as  column  0.  Elements 
k=141  and  168  were  the  two  elements  in  the  extreme  right  column, 
identified  as  column  13,  when  the  array  was  viewed  from  the  back.  The 
estimated  RF  phase  shifts  deviated  from  the  bench-measured  values  by  5*  or 
less.  The  average  deviation  was  slightly  less  than  2.5*.  This  agreement 
was  very  much  better  than  expected,  because  the  bench  measurements  were 
carried  out  with  a  differenct  CW  source  about  twenty  months  earlier. 

Since  then,  the  array  had  been  disassembled  and  reassembled  at  least 
twice . 


The  estimates  of  array-element  tuning  phases  were  examined 
indirectly  by  sweeping  a  beam  across  the  CW  source.  The  results  are  shown 
in  Fig.  7.  Here,  the  solid  curve  was  the  antenna  pattern  constructed  with 
estimated  values  of  8^$  and  The  broken  curve  was  constructed  with 

ideal  values  defined  as  ■  £x45*  and  $k  ■  0*.  There  is  no  significant 
difference  between  the  antenna  patterns. 

B.  Distorted  Antenna 

The  distorted  antenna  is  shown  in  Fig.  8.  The  feed  horn  was 
pointed  at  the  left  edge  of  the  array.  The  new  results  on  relative 
excitation  voltages  and  antenna  patterns  are  shown  in  Figs.  9  and  10, 
respectively.  In  Fig.  9,  the  excitation  voltages  on  the  left  of  the  array 
are  larger  than  those  on  the  right.  In  Fig.  10,  the  peak  of  the  antenna 
pattern  constructed  with  ideal  values  and  ^  is  on  the  left  of  the 
correct  position.  The  widths  of  the  peaks  in  the  two  antenna  patterns  are 
also  larger  than  those  in  Fig.  7. 


COLUMN  NUMBER 

Fig.  6  Estimates  of  relative  array-element  excitation  voltages.  The  dots 
denote  the  positions  of  the  defective  elements.  The  excitation 
voltages  at  these  defective  elements  were  less  than  102  of  the 
theoretical  values.  The  excitation  voltages  at  the  functioning 
elements  were  more  than  502  of  the  theoretical  values. 


Fig.  7  Antenna  patterns  constructed  with  estimated  (solid  curve)  and 

ideal  (broken  curve)  values  of  RF  phase  shifts  and  tuning  phases. 
The  ideal  values  are  defined  as  B^u,  *  lx45*and  ^  =  0*  for  all 
values  of  k  and  l.  The  antenna  is  distortion-free.  The  CW  source 
direction  is  marked  with  an  arrow  above  the  antenna  patterns. 


-15  D  15 
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Fig.  9  F.stimates  of  relative  array-element  excitation  voltages  in  the 

distorted  antenna.  The  dots  denote  the  positions  of  the  defective 
element  s . 
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10  Antenna  patterns  constructed  with  estimated  (solid  curve)  and 
ideal  (broken  curve)  values  of  RF  phase  shifts  and  tuning  phases. 
The  feed  horn  is  pointed  at  the  left  edge  of  the  array.  The  CW 
source  direction  is  marked  with  an  arrow  above  the  antenna 
patterns.  Note  that  the  peak  of  the  broken  curve  is  in  the  wrong 
posit  ion. 
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The  results  in  Figs.  9  and  10  were  consistent  with  the  distortion 
of  the  antenna.  The  feed  horn  was  pointed  at  the  left  edge  of  the  array. 
Therefore,  its  coupling  with  the  array  elements  on  the  left  was  stronger 
than  its  coupling  with  those  on  the  right.  This  explains  why  the 
excitation  voltages  on  the  left  were  larger  than  those  on  the  right.  The 
rate  of  decrease  in  excitation  voltage  from  left  to  right  was  rather 
rapid,  and  resulted  in  a  reduction  in  the  effective  aperture  of  the 
array.  This  aferture  reduction  was  the  main  reason  for  the  broader  beams 
shown  in  Fig.  10.  A  qualitative  explanation  for  the  position  of  the  main 
lobe  in  the  broken  curve  is  now  given.  In  Fig.  8,  the  opening  of  the  feed 
horn  was  farther  away  from  the  right  hand  side  of  the  array  than  the 
left.  The  compensation  for  this  longer  distance  was  essentially  the  same 
as  the  compensation  required  to  steer  a  beam  of  a  distortion-free  antenna 
to  the  left  of  the  boresight  direction.  Hence,  the  main  lobe  of  the  beam 
pattern  constructed  with  the  ideal  values  of  and  <t>^  was  on  the  left  of 
the  correct  position. 


VII  SUMMARY 

A  technique  to  measure  the  relative  array  element  excitation 
voltages,  tuning  phases,  and  RF  phase  shifts  of  a  phased-array  antennas 
has  been  described.  In  experiments  with  a  phased-array  antenna  at  the 
Communications  Research  Centre  in  Ottawa,  Canada,  the  estimates  of 
relative  excitation  voltages  were  found  to  be  in  good  agreement  with 
theoretical  values.  The  estimates  of  RF  phase  shifts  were  also  in  good 
agreement  with  bench-measu; ed  values.  Antenna  patterns  constructed  with 
the  estimates  of  RF  phase  shifts  and  tuning  phases  also  had  their  main 
lobe9  pointed  at  the  CW  test  source. 
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APPENDIX  A 


AN  ALGORITHM  TO  DETERMINE  THE  BEST-FIT  CIRCLE  THROUGH 
A  SET  OF  POINTS  IN  THE  ARGAND  DIAGRAM 

The  best-fit  circle  through  a  set  of  points  {Ak£;  1=0 , 1 ,2 , . . .  ,  7 } 
in  the  Argand  diagram  is  identified  here  as  the  circle  which  minimizes  the 
sum 


S 


i  (4 1  -  ak> 2 


1=0 


where  dk£  is  the  distance  of  A^£  from  the  centre  B^  of  the  circle  and  ak 
is  the  circle  radius.  Distance  dk£  is  given  by 

dk£=  (Ak£  "  Bk)2 

=  (x£  -  xb)2  +  ( y £  -  yb)2  »  (A. 2) 

where 

(x£,y£)  =  (Re (Ak£) ,  Im(Ak£))  ,  (A. 3) 

and 

A  A 

(xb»yb)  =  (Re(Bk) ,Im(Bk) )  .  (A. 4) 

Minimization  of  S  in  (A.l)  with  respect  to  ak  yields 

7 

“  Aak  l  [ ( x £  -  xb)2  +  (y£  -  yh>2  -  ak)  (A. 5) 

6ak  1=0 


=  0 


f 


which  leads  to 


where 


2xxb  +  2yyb  =  (x2  +  y2)  +  (xb  +  yb  -  ak) 


— nr 

x  y  = 


1 


(A. 6) 


(A. 7) 


Similarly,  minimization  with  respect  to  xh  and  using  (A. 6)  leads  to 


(A. 8) 


2x  2xc  +  2xyyb  =  (x3  +  xy2)  +  (xb  +  yb  “  a^x 
Minimization  with  respect  to  y»0  and  using  (A. 6)  leads  to 

2xy  xc  +  2y2yh=  (x2y  +  y3)  +  (xb  +  yb  -  a^)y  .  (A. 9) 

2  2  2 

Upon  elimination  of  (xh+yh-a^)  from  (A. 6),  (A. 8),  and  (A. 9), 

2(x2  -  xx)xb  +  2(xy  -  xy)yb  =  x3  -  xx2  +  xy2  -  xy2  ,  (A. 10) 

and 

2(xy  -  xy)xb  +  2(y  -  yy)yb  =xy-xy+y  -  yy  .  (A. 11) 

Solving  (A. 10)  and  (A. 11)  for  (x^,,^),  and  using  (A. 4)  leads  to  (16)  and 
and  (17)  Radius  a^,  calculated  with  (A. 6),  (16)  and  (17),  is  given  by 
(22). 

In  summary,  the  best-fit  circle  through  a  set  of  points  in  the  Argand 
diagram  is  identified  as  a  circle  which  minimizes  the  sum  S  given  by 
(A. 1 ) .  The  circle  centre  is  given  by  (16)  and  (17),  and  the  circle  radius 
is  given  by  (22). 
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APPENDIX  B 


COMPUTATION  OF  REFINED  ESTIMATES  {Yk£5  £-0,1 ,2 , . . . ,7} 


FROM  £=0,1, 2,..., 7} 


In  general,  the  improved  estimates  {y^IP;  4*0,1 ,2, .. .7}  given  by  (26) 
do  not  satisfy  (10).  This  appendix  is  concerned  with  the  construction  of 
refined  estimates,  denoted  by  {Yfc£»  £=0 , 1 ,2 , . . . ,7 }  which  satisfy  both  (9) 
and  ( 10) ,  i .e . , 


Yk£  "  ^k(£-4)  =  8k4,  £=4,5, 6, 7 


(B.l) 


Ykl  +  ^k2  =  Yk0  +  Yk3 


(B.2) 


Th^refined  estimates  are  linear  combination  of  the  Yk£  s.  If 
the  Yk£  s  satisfy  (10)  already,  one  has  Yk£  *  Yk£  for  all  values  of  £. 


This  property  suggests  as  trial  solutions 


Yk£  “  Yk£  +  e£6k,  £-0,1,2,.  ..  ,7 


(B.  3) 


where 


x  .  J2)  1(2)  J2)  J2)  x 

6k  ~  Ykl  +  Yk2  -  Yk0  ~  Yk3  +  Pk  x  360 


(B.4) 


and  e£  is  a  scaling  factor  to  be  determined.  Here  pk  is  an  integer 

which  restricts  6k  to  the  range  [-180*,  180*).  The  choice  of  this 
range  takes  into  account  that  <5k  is  zero  if  (10)  is  already  satisfied. 
From  the  symmetry  properties  of  the  subscripts  on  the  right-hand  side  of 
(B.4),  one  gets 


eg  =  -e  i  =  -e2  =  e3 


(B.5) 


Substitution  of  (B.3)  -  (B.5)  into  (B.2)  leads  to 


(l-4eo)<$k=0 


(B  .6) 


In  general  6^  *  0.  Therefore,  (l-4e0)  =  0,  giving 


e0  =  ~e i  =  -e 2  =  e  3  =  0.25  (B . 7 ) 

From  (B.l),  (B.3),  and  (26) 

A  A(2) 

e46k=  -  Yk£ 

/v  a  *  /  ry  \ 

=  fYk(£-4)  +  ek4^  ~  f  ^kC  £-4)  +  &k4l 
=  e4-46k  ,  (B  .8) 

£=4 ,5,6,7  . 

This  relation,  together  with  (B.7),  gives 

e4  *  -e 5  =  -eg  -  ey  -  0.25  .  (B.9) 

A 

In  summary,  improved  estimates  { Yk£ »  ^=0,1 ,2 , . . . ,7}  are  calculated 

a 

from  ;  4*0,1 ,2 , . . . ,7}  via  (B.3)  using  (B.4),  (B.7)  and  (B.9).  It  may 

be  noted  that  the  magnitude  of  &k  i-8  a  measure  on  the  accuracy  of  the  RF 
phase  shift  estimates.  If  j'Sknj  6*>  the  estimates  are  usually  very 
accurate.  However,  if  |6kj  90°,  the  errors  in  the  estimates  are  usually 
very  large. 


ARRAY  ELEMENT 


RF  PHASE  SHIFT  8k*  (DEGREES) 


k 

(ROW, 

COL.) 

ESTIMATED 

BENCH-MEASURED 

1=  1 

2 

D 

1 

2 

4 

121 

(  2 

,  o  ) 

47 

89 

166 

49 

94 

170 

148 

(  0 

,  o  ) 

49 

84 

175 

48 

88 

173 

175 

(-2  , 

0  ) 

48 

85 

167 

49 

90 

169 

141 

(  1 

,  13  ) 

48 

104 

189 

49 

104 

194 

168 

(-1  , 

13  ) 

45 

104 

196 

44 

103 

193 

TABLE  I:  Independent  RF  phase  shifts  obtained  by  estimation  and  by 
bench  measurements 
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